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Certain representations of the Cuntz relations, and a 
question on wavelets decompositions 

ov 

C^l . Palle E.T. Jorgensen 



Abstract. We compute the Coifman-Meyer-Wickerhauser measure for 
r \ , certain families of quadrature mirror filters (QMFs), and we establish that 

. ■ for a subclass of QMFs, contains a fractal scale. 

•5 

B: 

1. Introduction 

It is known that certain representations of the Cuntz relations (or equiv- 
alently the Cuntz algebras On [5] serve as a computational tool for wavelet 
analysis. The reason for this is that a representation of the Cuntz relations on a 
, Hilbert space 7i induces a successive subdivision of 7i into orthogonal subspaces. 

' To get the scaling feature of wavelets enter into our representations, build a fam- 

, ily of operators from the endomorphism z — > z N on the circle T, and N carefully 

chosen multiplication operators on L 2 (T) with respect to Haar measure on T. In 
the dyadic case, we work with z — ► z 2 and with two multiplication operators, 
yS^, ' i.e., two functions on T. In the language of signal processing, these functions 

, m , and mi are called low-pass and high-pass filters, respectively; or the pair 

is called a quadrature- mirror system of filters. By a standard procedure [7], we 
construct the scaling function (or father function) tp on K from vtlq, and the 
^ ' wavelet function (or mother function) ip from mi. We get two corresponding 

isometries which yield a representation of the Cuntz algebra Oi- These repre- 
' sentations have been studied in earlier papers, see e.g., [T | |2[ l3 l l4l 15 1 112] . The 

problem in wavelet theory is to build orthonormal bases in L 2 (R) from these 
data. This can be done [7], and the wavelet bases have advantages over the 
earlier known basis constructions, one of the main advantages is efficiency of 
computation. 
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While choice pyramids (see e.g., |15) and |16|) have been much used in 
computations already, the presence of representations of On yields a Hilbert 
space geometric way of realize such combinatorial pyramids. Representing the 
integers No = {0,1,2,...} using the Euclidean algorithm, and using a corre- 
sponding pyramid, we get a orthogonal basis for L 2 (R), ipo, ipi, ip2, ■ ■ ■ called a 
wavelet packet. A wavelet packet decomposition arises by a combinatorial selec- 
tion of a subset of the index set No for these basis function, and a corresponding 
subset of the scaling operations x — ► N^x (i.e., a selection of a subset of the 
integers Z) such that this selection of subsets yields an orthonormal basis for 
L 2 (R). As pointed out in |15) and |16| . the selection of such feasible pairs of 
subsets involves a combinatorial tiling question, and a problem in measure the- 
ory. In this paper we show how these questions may conveniently be addressed 
with the use of the non-commutative harmonic analysis coming from the theory 
of representations of the C*-algebra On- 

And in the course of the paper, we show that the scope of our method is 
wider than the original context of standard wavelets: It applies equally well to 
the class of fractals that arise from affine iterated function systems (IFS); see 
e.g., [E], US], and [Hi- 
Wavelets and wavelet packets are orthonormal bases, or Parseval frames, for 
the Hilbert space L 2 (R) which are built from a small set of generating functions 
by applying only two operations to them, translation by the integers Z C R, 
and by dyadic scaling t — ► 2H,j 6 Z. The best known case is that of singly 
generated wavelets. Then there is a function ijj £ L 2 (R) such that the family 

(1.1) l2 j / 2 ib(2H-k)] 

1 J i.fcez 

is an orthonormal basis for L 2 (R). Wavelet packets (in the orthonormal case) 
are sequences (yn) ng No an( ^ subsets J c No x Z such that the two families 

(1.2) {2 J/ V (2H - k) | (n,j) eJ,kez] 
and 

(1.3) {(p n (t-k) | neN ,/ce Z} 

are both both orthonormal bases. 

In the seminal paper [5], the authors Coifman, Meyer, and Wickerhauser 
proposed a selection of systems l|1.2|l from 1)1.30 : the idea being that 1)1.3(1 may 
be constructed by an effective and relatively simple algorithm. The construction 
starts with a quadrature-mirror wavelet filter, i.e., a Fourier series 

(1.4) m (x) = Y,a k e- i2 * kx 

fcez 

such that 

(1.5) ^2 ak a k+2e = Sa,e 
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and 

(1.6) ^a fe = V2. 

If the expansion in ljl.4jl is a finite sum, then the functions ipo, ipi, tf2, • ■ • will be 
of compact support. The advantage with using (|1.2(l over l jl.lt or Ijl.Hjl is that 
the basis functions in (jl.2jl are better localized in time-frequency, as is spelled 
out in [5] and |15| 

The authors of |5] propose a certain integral decomposition of L 2 (R) which 
is based on a certain measure /io on the unit-interval [0,1]. The idea is the 
following: The quadrature rules (jl.5jl lead to a dyadic decomposition of the 
Hilbert space £ 2 (Z), repository for the wavelet coefficients: For each k £ N, and 
each dyadic interval 

(i-7) [^+2- fc M=|+---+!, 

we assign a closed subspace H (£) c £ 2 , such that the 2 k distinct subspaces H (£) 
are mutually orthogonal, and 

(i.8) y, h & = £2 

If eo denotes the vector eo (j)'= #o,j i n an d f(£) the orthogonal projection 
onto H. (£) , then the measure /j,q is determined by 

(1-9) Al0 ([e,e + 2- fc )) = ||P(e)e || 2 

where ^ is the dyadic rational ljl.7jl . 

More generally, if / € H, ||/|| 2 = 1, we set fif (•) = ||-P(-)/l| 2 ; an d note 
that each fif is a probability measure. These measures are shown in section [3] 
to dictate the wavepacket analysis in the Hilbert space L 2 (R). 

The expectation is that /io is absolutely continuous with respect to the 
Lebesgue measure on [0,1]. We give a formula for /i directly in terms of the 
coefficients (aj) eZ in ljl.4jl . and we give a one-parameter family of coefficients 
cij (/?) where [3 parameterizes the circle, and where for \{3\ < 7r/4, the measure 
Ho on "small" dyadic intervals J is governed by the formula 

(1.10) (dyadic interval J) ~ | J\ s , 

where |J| denotes the length of J, and where 

In' — 2 



(1.11) s = - V 



In 2 

In the interval for (3, the coefficients aj (f3) are real-valued, and when \/3\ < 7r/4, 
a 09) > 1/V2- 
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2. Fourier Polynomials 

Let T = {z G C | \z\ = 1}, and let D € N. Consider the following two 
functions 

2D-1 

(2.1) m 0) = a k zk 

k=0 

and 

_2D-1 



fc - fe 
0-2D-l~kZ 



(2.2) mi(2) = ^ -1 moH 

2Z»— 1 

= E (-D 

fc=0 

We shall need several lemmas. 

Lemma 2.1. The finite sequence ao, a±, ■ ■ ■ , 02/3-1 in ( 12.1)) satisfies 

(2.3) ^ afeafe+2£ = £0/ 

if and only if anyone of the following equivalent conditions holds: 

(a) The matrix function 

D-l 

(2.4) U(z)=J2A k z k ,zeT 

k=0 

is unitary where the coefficient matrices are 

(2.5) A k =( a2k a2k+1 ) forfe = 0,l,--- ,D-1. 

V a 2(D-fc)-l ~ a 2(D-fc-l) / 

(b) The operators 

(2.6) (Sif) (z) = ra, (2) / (z 2 ) ,/eI 2 (T),i= 0, 1, 
and their adjoints S 1 * satisfy 

ic, 7) / STSj = kjl, i,j = 0,1 

(c) The following matrix function 

(2.8) M (z); = ^f mo ^ mo 5" Z ! VzeT, 

y/2 \ mi (2) mi (-z) y 

is unitary i.e., M (z)* M (z) = I, for all z 6 T. 

Remark 2.2. In the summation (|2.3|l . it is understood that the terms are 
zero if the summation index is not in the range which is specified for the sequence: 
If for example D = 2, then the conditions (|2.3)l spell out as follows: 



(2.9) 



|ao| 2 + |ai| 2 + M 2 + |a 3 | 2 = l 
a a 2 + aia 3 = 
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The interpretation of (|2.4|) is that the matrices Aq, Ax, ■ ■ ■ Au—x are Fourier 
coefficients for a matrix function z — > U (z) defined on T. The requirements is 
that for each z £ T, the 2 by 2 matrix U (z) is unitary, i.e., that 

(2.10) U{z)*U{z) = I, 261 

In general, for a function of the form l|2.4|l . the unitarity condition (|2.10|) is 
equivalent to the following conditions on the constant matrices Aq, A\, ■ ■ ■ : 

(2.11) J2 A 

fe 

or equivalently 

(2.12) J2 A k A k+n=6o,n- 

k 

Again if D = 2, these conditions flesh out as follows: The two matrices 

(2.13) A 
satisfy 

(9 xa) / A a AZ + A l A\ = I 

^• MJ \ A X A* Q =Q 

Proof, of Lemma 12.11 The coefficients ao,ai, • • • ,a 2 D-i ar e given and 
the two functions m (z) and mi(z) are defined as in (|2.1|) and i|2.2[l . Even 
though they are defined for all z 6 C, we shall use only the restrictions to 
T. It is convenient in some later calculations to introduce the substitution z\— 
e -i27T£c. an j w j^ n this substitution we shall also view mo and ra\ as one-periodic 
functions on the real line R, and make the identification T = K/Z. We now turn 
to the implications: (|2.3|) =>fa N ). We already noted that the unitarity of the 
function U(z) in (|2.4|) may be expressed by the conditions (|2.11|) . When l|2.5|l 
is substituted in l|2.11|l . it is immediate that the two identities l|2.3|l and (|2.11() 
are equivalent. (a)=>(b). Using the normalized Haar measure /i on T and the 
corresponding inner product 



( a 




\ and Ax = I 




03 \ 


V a 3 


-02 j 




v ax 


-ao J 



(2.15) (f\g)--= / f(z)g(z) dn(z) 

Jt 

of the Hilbert space L 2 (T), we get the formula 
(2.16) 

(S*f) (z) = - m i H/ ( w ) ' dcnncd for / e L 2 (T) , z e T , and i = 0, 1. 

As a result, we get the formula 

(2.17) {S*Sjf)(z) = ± ]T mjwjmj (w) f (z) ; 
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which states that the four operators S*Sj are multiplication operators. Hence 
the first part of 1|2.7[1 reads: 

(2.18) ~ ^2 m i ( w ) m j ( w ) = , z e T. 

w£T,w 2 =z 

Substitution of l|2.1[l into (|2.18|) shows that these conditions amount to the same 
sum rules (|2.3|l . or equivalently (|2.11() . The second equation from 1|2.7[) may be 

restated as 

(2.19) \\S* f\\ 2 + \\S*J\\ 2 = Il/H 2 , for all f e L 2 (T) ; 
or equivalently 



l 



(2.20) 



- V m l (w)f(w) d/j,(z)= / |/(z)| 2 dfx(z). 
z , jt 



2 

—z 



i=0 

But the summation on the left-hand sides fleshes out as follows: 
l 



\1>2 E ^Hm^o/H/CO 

i=0 w 2 =Q 2 =z 

= \ E UE^H^(C))/H/(C) 
= ^ E <W/M/(0 



2 - ■'"•' ,2 



and 

(2.21) /i ^ |/H| 2 ^(z) = / |/(z)| 2 d M (*) 

w 2 —z 

by an elementary property of the Haar measure. This proves (a)=^(b); and in 
fact equivalence. (b)=4>(c). This implication also goes in both directions, and 
it is implicit in the previous step where we showed that the identities (|2.18ll are 
equivalent to l|2.3fl . But an inspection of matrix entries shows that (|2.18l) is a 
restatement of the unitary property for the matrix M (z) in (c). □ 

Corollary 2.3. Consider a Fourier polynomial 

2D-1 

(2.22) m{z) = a k zk 

k=0 

as in Lemma \'2. II and let 5* = S m be the operator on L 2 (T) defined by 

(2.23) (Sf)(z)=m(z)f(z 2 ) , z E T. 

Then the following conditions are equivalent: 
(i) S is an isometry; 
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(ii) \m (z)\ 2 + |m (— z)\ 2 — 2 , for z G T; and 

(iii) J2k®kak+2e = So,e- 

Proof. The equivalence (i)<^=>(ii) is immediate from an inspection of (|2.18|) 
in the special case i = j = 0. To see (ii)<^=S>(iii), it is simplest to write out the 
Fourier coefficients of the function 



(2.24) z^ 1 - Yl I^HI 2 - 



A substitution of liFT^, into $171^ yields 



?EE E w" 



J K W 2 =Z 



E E {\ E ^ 

n j,k s.t. k—j—2n \ w 2 =z J 



n j,k s.t. k—j—2n 

= E E "j a j+ 2 ™ 



and it follows that the sum in 1|2.24(1 is the constant function 1 if and only if 
(iii) is satisfied. Note that (iii) is a restatement of condition (|2.3|) in Lemma 
In the following we consider operators S on the form (Sf)(z) — m(z)/(z 2 ), 
z € T, / e L 2 (T), and their adjoints. But we will need these operators realized 
explicitly in the sequence space I 2 = L 2 (T) . □ 



Lemma 2.4. . 

(a) The matrix representation for the operator S in Corollary 12.31 relative to the 
standard basis e„ (z):= z n , n € Z , for L 2 (T), is 

(2.25) (S£) n =E a «-2^/ 



and the adjoint S* is 

(2.26) (S*0n = E S *-a»&- 
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(b) The 2D-dimensional subspace C spanned by {efc | —2D + 1 < fc < 0} is in- 
variant under S*, and the corresponding 2D by 2D matrix is 



(2.27) 



«2 
&2D-2 





"3 




«o 

«2 






















02D-1 


«2D-2 


a 




















0,2 


«i 


no 













0,3 


0,2 


a\ 








0>2D-2 




















a-2D-\ 


0-2D-2 


02D- 




















0.2D- 



V 

Specifically, the case D = 2 is 
(2.28) 



i / 



/ ao 

0,2 Oi 





\ 







a 

a 3 a,2 ai 
a 3 / 



and the space C is then spanned by {l, 
(c) For all neZ, there is a k £ N such that 



— 1 —9 

z ,z ,z 



(2.29) 



5* fc e„ e C. 



Proof. The matrix representations 12.25fl . and l|2.26|) . follow directly from 
the formulas (|2~23|) and (|2~Tol) . To prove (b), set 

J (D) : = {fc S Z | -2D + 1 < fc < 0} , 

so that 

£ = span {efc | fc e J (D)} . 
We claim the following implication: 

(2.30) & = for fc e Z\ J (D) =>• ^ a fc _ 2 n& = for n G Z\ J (D) . 

fc 

The invariance of C under S* is immediate from this. Now suppose = for 
fc G Z\J(D). Then one of the factors in afc_2n£fc vanishes if n G Z\J (D); and 
the implication (|2.3U|) follows. To understand (c), iterate the matrix formula 
We get 



(2.31) 



S* k e r . 
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where the summation on the right-hand side is overp S 1, i\, ■ ■ ■ G [0, 1, ■ • • , 2D- 
1], subject to 

h + 2i 2 + h 2 fc ~ 1 i fe + p2 fe = n, 

so the range for the p-index is roughly divided by 2 with each iteration of S* on 
the basis vector e„ (z) = z n . □ 

Remark 2.5. Conclusions (a)-(b) in Lemma f2.4l hold for the smaller sub- 
space M.:= spanjefc | —2D — 2 < k < 0}. The matrix of the restricted operator 
5*1^ is then obtained from i|2.27[) by deletion of the last row and last column. 
For example, the reduced matrix corresponding to l|2.28ll is 

/ a 
F: = I a,2 ai ao 
\ a 3 a 2 

However, property (c) in Lemma [2.41 is not satisfied for M.. To see this, note 
that D = 2, and C = M. ® Ce_3 in the example. As a result 

(2.32) ^e_ 3 | S*V 3 J} = (a 3 ) k for all k e N. 

If 03 7^ 0, it follows that S* e_3 is not in Ai, no matter how large k is. 
3. Measures Induced by Representations of the Cuntz Relations 

The system of coefficients studied in section [21 give rise to representations 
of the Cuntz relations. We recall the construction briefly; see [5] for more 
details. When the coefficients a , ai, • • • are given, satisfying l|2.3|l . then the two 
operators So and Si from Lemma l2~T1 fb) satisfy the relations 1)2. 7f) . We say that 
the two operators satisfy the Cuntz relations, or equivalently that they define a 
representation of the Cuntz algebra O2 on the Hilbert space H — L 2 (T) . The 
norm-closed algebra of operators on Ji generated by So, Si, Sq, and SI is known, 
[S] to be a simple C*-algebra; and it is unique up to isomorphism of C*-algebras. 

We will also need to Hilbert space L 2 (M.) . If in addition to (|2.3(l , the numbers 
(dfc) satisfy 

2D-1 

(3.1) Yl ak = A 

fe=0 

then it is known [5] that the following system of three equations has solutions 
inL 2 (K): 

(3.2) (p 2n (x) — V2 a k <pn (2x - fc) , 



(3.3) (fi2n+i (x) = V2^ (- 1 ) fe a2D-l-kfn (2x - k) 

k 

and 

(3.4) / tpo (x) dx = 1. 
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Moreover, for every / G L 2 (R), the following Parseval identity holds: 

OC r. 

(3-5) EEK^(-- fc )l/)l 2 = ll/H 2 = / l/^)! 2 dx - 

When a system, such as 

(3.6) Wn(--k) |neN ,fceZ}, 

satisfies (|3.5f) we way it is a Parseval frame, or a normalized tight frame. There 
are simple further conditions on the coefficients (afe) each of which guarantees 
that (|3.(i|) is in fact an orthonormal basis (ONB) in L? (R), i.e., that 

(3.7) (tpn (• - k) | (fi m (• - I)) = S n , m Sk,e- 

To understand how to select subsets J C No x Z such that the set of functions 

(3.8) F(J): = fcltpnfix-k) | (n,j) e J, fc e Z j 

forms a Perseval frame, or an ONB, for L 2 (R), the authors of [5] and |16| sug- 
gested a family of measures associated with the operator systems <|2.7ll . We 
proved more generally, in |14j that these measures are associated with any rep- 
resentation of one of the Cuntz algebras On, N = 2, 3, • • • . 

The idea is simple: Let (<Si) i=0 be two operators in a Hilbert space TL such 

that 

l 

(3.9) S*Sj = S id I and E = I 

i=0 

hold. Then for every multi-index 

(3.10) £(*i,*2,--- ,t fc ),*j G {0,1} 
the operator 

(3.11) P^):^---^* •••£* 
satisfies 



p(0* = -p(0 = -P(0 



2 



Moreover 



(3-12) E P ^) =/ 
and 

(3.13) P(0P(v) = S tn P(0 

where the summation in (|3.12l) is over all multi-indices of length k, and in (|3.13() 
the multi-indices £ and r\ have the same length. We say that, for each k £ N 
we have a dyadic partition of the Hilbert space TL into orthogonal "frequency 
bands." 
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Theorem 3.1. Let m .'T — > C be a function which satisfies m (1) = \/2 

and 

(3.14) |TO (z)| 2 + |m (-z)| 2 = 2 forzeT. 
Set 



(3.15) mi (z) = z mo (— z), for z £ T; 

and let Si, i = 0, 1, be the operators 1|2.6[1 . Suppose the sequence tpo, cpi, cp2, ■ ■ • 
defined by 

(3-16) <p^ (x) = -^ro (|) fa (|) 

(3-17) <Pito+i(x) = -^=rni (|) fa (|) 

for n = 0, 1, • • • , x e R defines an orthonormal basis in L 2 (R), where we set 

m 3 (x):= m 3 (e l27ra; ), j = 0, 1. Then 

(3.18) 

i i 

2%ip„ (2 p t - k) = ^ ■■ ■ ( e i I S i P " ' S h e k) V2Pn+ii+i 2 2+-+i p 2P-i (t - j) 

i 1= Q i p= 0j£Z 

holds for all p, n £ No and k £ Z. Moreover a subset J C No x No has the 
corresponding set 

(3.19) {2« ip n (2 p t - k) | (p, n) £ J, k £ z} 

define an orthonormal basis ( ONB ) in L 2 (R) if and only if the sets [2 p n, 2 P (n + 1)) 
form a non-overlapping partition of No- 

PROOF. Since the functions in (|3.16(l - (|3.17|) form an orthonormal basis for 
L 2 (R), we have 

(3.20) Y fa (x + j)fa7 (x + j) = 8 n y for n, n' £ N , 
and 1|3.18[) will follow if we check that 

(3.21) f tp m (t-j)2Sip n (2 p t -k) dt= ( ej \ ST ■ ■ ■ S* e fc ) 

if m = 2 p n + i\ + i%2 + ■ ■ • + i p 2 p ~ 1 , and zero otherwise. But H3.21[l follows 
from Fourier duality, and a substitution of formulas 13. 16(1 - 1)3. 17|l and l(3.2()|l . 
We leave the computation to the reader. The second conclusion regarding the 
basis properties (ONB) of the functions in 13.19fl follows from (I3.18|l when we 
note that 

m = 2 p n + it + i 2 2 H h i p 2 p ~ 1 

ranges over [2 p n, 2 P (n + 1)) as i\, ■ ■ ■ , i p vary over Z 2 x • • • x Z 2 . □ 

p times 
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Example 3.2. The best known special case of (|3.19|) is when 

J = {(0,1), (1,1), (2,1), (3, l),---}, 

and the corresponding partition of No is simply [2 P , 2 P+1 ) for p = 0, 1, 2, • • • . In 
this case, we only need the first two functions from the sequence (p n of i|3.16|) - 
(13.17(1 . i.e., (fo — ip, and <p\ = ip. The other extreme is J = {(0, 1), (0, 2), (0, 3), 
• • • }. In this case, each section in the corresponding partition of No is a singleton, 
and all the functions {ip n | n E No} from ((3. 16(1 - ((3. 17(1 are needed in the ((3.19(1 
family. However, no dyadic scales are used. Between these two extremes there 
is a countable infinite family of non-overlapping partitions corresponding to the 
choice specified in ((3. 19(1 : for example, the pairs (p,n) in some J from ((3.19(1 
may be specified as in the following table (Table 1). 



Table 1. 

A set J consisting of points (p, n) as follows: 
p2222444666888 
^0123123123123 



In this case, we are using only the four functions ipo> </?i> </?2, and tp 3 from ((3.1611 
(13.17(1 . and the segments from 1(3.19(1 are {0, 1, 2, 3} followed by 

[2 2fc j, 2 2fc (j + 1)) 

where k — 1,2,3, and j £ {1,2,3}. As a result, all the dyadic scaling 
numbers 2 2fc are used. The resulting ONB in 1(3.19(1 consists of the following 
functions, 2ip Q (2 2 t - £ ) , 2 k ipj (2 2k t - t) , where £ ,£ € Z, k = 1, 2, • • • , and 
J e {1,2,3}. 

Since {ip m (• - j) | m € N , j G 1} is an ONB for L? (K), every } E L 2 (M), 
|| /|| 2 = 1, defines a probability distribution Pf on No X Z by 

(3-22) Pf(mJ): = \(<p m (--j)\f)\ 2 . 

Corollary 3.3. Let p,n e No, k € Z, and let ^ 2P {.;•) be the 

corresponding probability distribution from 1(3. 22() above. Then the marginal 
distribution in the first variable is 

(3.23) S m 2Pn+i 1 +i 2 2-\ h« P 2P- 1 ^e k 



2 2p ' 2 



+ 



2i' 



1 

2p 
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PROOF. In view of 13.2211 and Theorem 13. II the marginal distribution is 



22 tp„(2P—k) 



(™>,j) 



Sm,2Pn+i 1 +i 2 2+---+i p 2P- 1 | \ e J I ^ip ' ' ' ^n ek 

2 



,2Pn+ii+i 2 2H \-ip2P~ 1 



= ^m,2"n+ii+i 2 2H hi P 2P- 1 Met 

which is the desired conclusion. 

The next result follows from or from |14| 



2?' ~2 



_p_ I __ 

2 p 2P 



□ 



Lemma 3.4. There is a unique Borel measure E defined on the unit interval 
[0, 1] such that 



(3.24) E 



2 k ' 2 



2 fe 2 k 

Moreover the measure E = Ep takes values in orthogonal projections in 7i, and 



(3.25) 



/ E (dx) = I 
Jo 



(3.26) 



E (Bi) E (B 2 ) = if Bi are Borel sets with B x Q B 2 



Remark 3.5. The expectation is that the measure class of E is that of the 
Lebesgue measure if E is derived from a wavelet representation. For the wavelet 
representations, the Hilbert space H is L 2 (T); and for a generic set of wavelet 
representations eo is a cyclic vector; see |13j . Hence, to check if E is in the 
Lebesgue class, it is enough if we verify that the scalar valued, positive measure 



(3.27) 



Mo (•): = (e \E(.)e ) = \\E (.)e \\' 



is absolutely continuous with respect to Lebesgue measure on [0,1]. Using 13.2411 . 
we see that the measure /io is given on dyadic intervals by the formula 



(3.28) 



+ 



i_k_ 

2 fc 



1 

2^ 



\Sie \ 



Part of the reason for expecting that the measure A*o (') = \\E (0 e ol should 
be absolutely continuous with respect to Lebesgue measure on [0, 1] rests on 
the observation that A = is always in the spectrum of the matrix l|2.27|l . If 
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A = is also a dominant eigenvalue, in a sense which we make precise below, 
then the asymptotic formula 



(3.29) 



u iu ii 1 1 i s 



would follow. From this the absolute continuity follows as well. But a closer 
examination of the matrices Fi,i = 0,1, shows that A = is not always a 

dominant point in the spectrum of Fq = Sq L, . 
We now show that is in the spectrum. 

Lemma 3.6. Let the numbers do, a%, ■ ■ ■ , &2_d-i satisfy (|2.3() and i|3.1[l . and 
let So and Sq be the corresponding operators on £ 2 ; see (|2.25|l and (|2.26|l . Then 
A = i is an eigenvalue for the finite-dimensional operator F = Sq \ m . If 
Fq is represented as a 2D — 1 by 2D — 1 matrix as in 

then the row- 
vector w — (1, 1, • ■ ■ , 1) is a left-eigenvector for Fq, or equivalently w* is a 



2D-1 



right-eigenvector, with the same eigenvalue for the adjoint matrix Fq*. 

Proof. Substitution of too (1) = y/2 into |to (z)\ 2 + \m (— z)\ 2 = 2 from 
Corollary 12.31 yields too (—1) = 0. Hence, z + 1 divides the polynomial too (z). 
If 

to (2) = (z + 1) C (z) , and 

2D-2 

C (z) — £iZ l , then it follows that 



i=0 
2D-2 



r( i) = £ A = v* * 
^ 2 V2 



As a result we get 

!a = £q 
a,i=£i + h-\ for 1 < i < 2D - 2 
02D-1 = ^2D-2- 

It follows from this that 



2 = 



and 



D-l j 



V2" 

Working out the matrix product, this is a restatement of wFq — -^w, or equiv- 
alently Fq*w* = iw*. □ 
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Theorem 3.7. Let a 0l ai, ■ ■ -a.2D-\ be as stated in Lemma 13.61 suppose 
|a | > max{|A| | A e spec (F ) \ {a }} , 

and set 

_ Ha | 2 
S ~ In 2 ' 

Then every non-empty open subset V in [0, 1] contains an infinite sequence of 
dyadic intervals J such that 

< hminf ( < hmsup ( ^fi) < oo. 
.rev V \J\ J ~ Jcv V \J\ J 

We say that the measure /Iq contains s as a fractal scale. 



PROOF. The details follow by combining Lemma 13.61 with the technical 
lemma which we state and prove in Section [5] below. □ 

Proposition 3.8. Let the numbers oo, • • • , £t2£>-i satisfy the two conditions 
(|2.3(l and H3.1(l and let So, Si be the corresponding two operators; see l|2.6|l - H2.7f) . 
Let /i (■):= \\E (•) e j| 2 be the measure (13.281) . Then 
(3.31) 



/'-o 



h , . h h . . h . 1 \ \ ^ i ,2-#(i=o) i i2-#(i=i) 



Remark 3.9. We will see later how this estimate yields information about 
a fractal component of /xq in the special case D = 2. 



Proof. Let F,:= S* 



, i = 0, 1. For k € N, and 7 ik € {0, 1}, set 

M 



£ = y H h ^ and F^:= F ifc • • • F h . Then F *e = a e and F-feo = a 2 D_ie . 

To see this, recall that 



2D-1 



and 



mi (z) = a 2 D-i - a 2 D-2Z H a z 



A10 ([^e + 2-' £ )) = ||^ e o|| 2 >|( e o|^e )| 2 
= l(^---^>o|eo)| 2 
= | a o| 2 - #(l=0) |a 2D ^| 2 - #(i=1) 



□ 



Lemma 3.10. (The case D = 2.) 
(a) For the case D = 2, the real valued solutions ao, ai, 02,03 to the three 
equations 

a o + a i + a 2 + a i = 1 
(3.32) <( a a 2 + aia 3 = 

a + ai + a 2 + a 3 = y/2 
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arc 



(3.33) 



a ° = 273 (1 + V2cos/3) 

01 = 2^3 (1 + V2sin/3) 

02 = ^(1-^2 008/3) 

a 3 = ^(l-x/2sin/?) 



for /3 e K. 

(b) The spectrum of the matrix 
(3.34) F = 



a 
a 2 ai a 
a 3 a 2 



obtained by restriction of S* to M — spanjeo, e_i, e_2} is 
(3.35) f~ 1 ™P -<*»(>' 



(c) The following two conditions are equivalent 

1 



(3.36) 

and 
(3.37) 



a > —= > max {| A| | A e spec (F) \ {a }} 
v2 



I"! < T 



(d) The eigenspace for the point a in the spectrum of F is spanned by the 
vector (*) where v £ R 2 is the solution to 

-l 



(3.38) 



ai a 
a.3 a 2 



«2 





(e) For the case of Daubechies's wavelet, we have 

a 



(3.39) 



(3.40) 

and 
(3.41) 

In particular, 
(3.42) 



ai = 



1+V3 
_ 3-jVj 



a 2 



a 3 = 



[3 = arccos 



iV2 
3-V3 

_ 1-V3 



iV2 



2V2 



slH ' r(F,= ' ao '71'271 



a < —= = max{|A| | A <E spec (F)} . 
V2 
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PROOF. Part (a) follows directly from section 4 in 2 . To compute the 
spectrum of the matrix F, note that its characteristic polynomial is 

(3.43) (A - a ) (A 2 - (ai + a 2 ) A + a\a 2 - a a 3 ) ; 

and when (|3.33|) are substituted, we see that the roots are as listed in l|3.35|l . 
The conclusions listed in parts (c) and (d) are immediate consequences of (b). 
Finally (|3.38|l follows from the observation that F has the form 



a 
OS.4-1) F= I a 2 G 

with G = ( ai a ° | and the fact that an is not in the spectrum of G. The 
\ a 3 a 2 J 

conclusions listed in (e) for Daubechies's wavelet follow by an application of 
(a)-(d). □ 

To better appreciate the geometry of the formulas (|3.33() . the reader may 
find Fig. 1 useful. 



The coefficients a (/3) , a x {(3) , a 2 (/?) , a 3 (/3) of (pO^I 

Figure 1. 

The four coefficients a^,a\,a 2 , a 3 are described by the two circles in Fig. 1, 
viz., by 
(3.45) 

l A 2 ( l A 2 l / l A 2 / l A 2 l 

a - — = +[a 3 - — = =7, and by ai - — -= + a 2 



2V2J V 2V2J 4' V 2V2J V 2V2J 4 ' 

This representation also makes it clear that, if one of the four coefficients van- 
ishes, then so does a second from the remaining coefficients. The resulting four 
degenerate cases are as sketched in Table 2. 
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Table 2. 

The degenerate cases when one (and therefore 
two) of the coefficients ao, ax, 02, 03 vanish: 



1+2 


1+z 3 


z+z 2 






V2 


V2 


V2 


z 2 -z 3 


1-z 3 


-z+z 2 


1-;. 


V2 




V2 


V2 


7T 

4 


7T 

4 


3tt 
4 


3tt 
4 



There are four distinct configurations, and all four are variations of the Haar 
wavelet, see Fig. 2. The second column in Table 2 is called the stretched Haar 
wavelet; see [3]. For all three, we have the formula 



(3.46) 



/'n 



h ik h ik 1 

~2 + " ' + ¥' ~2 + " ' + ¥ + ¥ 

which proves that fj,Q is the Lebesgue measure on the unit-interval [0, 1] for all 
four variations of the Haar wavelet. The two matrices Fq and F% corresponding 
to the four variants of the Haar wavelet filters in Table 2 are as follows: 



1 



10 

1 Its 
010 



000 
1 1 



000 
1 



)*, 
001/ v 1 1 




1 







Its 
0-11 



1 

1 )73' 
0-10' 





1 







1 )*, 

001/ V 



1 
-1 1 
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The father function ip and the mother function %p for the Haar Wavelet 
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Figure 2. 

Remark 3.11. We shall need Lemma f3. 61 for the two matrices 

a \ / a 3 

(3.47) F = | a 2 a\ a and Fj, = a% -a 2 a 3 
a 3 a 2 / \ -ao a3 

For the range \j3\ < ?, we have <|3.36[1 satisfied by spec (i^o)- But then the point 
fl 3 = a 3 (P) satisfies | (Z3 j < ;4=, and 

(3.48) —= = max{|A| | A e spec (Fi)} . 
V2 

Theorem 3.12. Let f3 G R, |/3| < ~, be given. Let a , a l7 a 2 , a 3 , be the 
corresponding numbers in (|3.33|l . Let (Si)l =0 be the operators in l|2.6(l . and 



(3.49) Fj: = 5* 

the matrices (I3.47|) . Let 



M 



(3-50) i = 77 ■ 



— H h — 

2 + + 2 h 



be a dyadic fraction. Then 

(3.51) lim a^ n F^F lk ---F n e = a* {l=0) af {t=1) (e + v) 
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where v is the vector in 
such that 



(3.52) 



Mo 



Moreover, if e > is given, there is an no £ 



,2n„2-#(i=0) 2-#(i=l) 



for all n > uq. More generally, given e, < e < 1, there is no € N such that 
#(*=o)^#(*=i)^ 2 < Mo 2^nr)) 



v\ — e 



(4 



Jin 



for all n 



< (4 (i = 0) fl 3 #(i=1) ) (i + IHlN 

Proof. Recall that eo + v is the eigenvector satisfying 

(3.53) F (eo + u) = a (e + u) , and (e | e + v) = 1. 
Moreover, we have 

(3.54) i^ *e = a e and F*e = a 3 e . 
Using l|3.36|) for Fq, we get 

(3.55) lira cu n F?w = (e | w) (e + v) 

n >oo 

for all w £ M. Applying this, and H3.54JI . to w:= Fi k ■ ■ ■ F^eo, the conclusion 
(TT5T|) follows. Since 



(3.56) 



Mo 



<5 o^eo 



(^••■^eolf , 



the second conclusion l|3.52ll follows from the first. 



□ 



Corollary 3.13. Let \(3\ < j, and let clq, a\, a 2 , 03 be the corresponding 

numbers in l|3.33[) . If the measure Mo (•) = ll^(') e o|| 2 is absolutely continu- 
ous with respect to Lebesgue measure dx on [0, 1]. Then the Radon-Nikodym 
derivative dfio/dx = / is unbounded in every open (non-empty) subset of [0,1]. 

PROOF. Let V C [0, 1] be a non-empty open subset. Then pick £ as in 
5D|), and n\ £ N such that 



(3.57) 



Then pick e, < e < |, and no > i%i such that (|3.52() holds for all n > uq. If / 
were bounded in V with upper bound C; then (|3.52l) implies the estimate 



(3.58) C-2~ n >al n a 2 * {l = 0) a 2 3 * (l=1) 



11 2 



for all n > tiq. Since a§ > \ , lim 



an ~2n 
n >oc * w 



co, which contradicts l|3.58|l . □ 
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Example 3.14. We introduced the constants a , aj., a%, 0,3 for Daubechies's 
wavelet in part (e) of Lemma 131)1 An inspection of l|3.33[l and (|3.39|) shows that 

f3 = arccos {j^^j ■ Hence 

sin/3 — cos/3 1 



2 2y^' 
and the three numbers in the spectrum of Fq are as follows: 

(3.59) ^<"<^ 
One checks that the vector 

(3.60) v*:=(l + V3)e-i + (l-V3)e- 3 



spans the ^-eigenspace for Fq. Since ao < the eigenspace 

1 

7=2 



(3.61) (weM\ F w = -^=w 



is orthogonal to eo- 

Corollary 3.15. Let /3 — arccos y ) 1 so * na ^ the numbers do, ai, 0,2, 
03 define the Daubechies wavelet. Let v* be the eigenvector in l|3.60|l . and pick 
w G M such that 

(3.62) F*w = -iffl, and (w | v*) = V8. 

v2 

Let £ be a dyadic rational given as in l|3.50|l . Then, for every e, < e < 1, there 
is no € N such that 

(3.63) I («, I F ik ... Fil eo) | 2 - e < MMj^Elj 

< I (to I F ifc • • • Fij e ) | 2 + e for all n>n . 

Proof. The argument follows the proof of Theorem 13.121 The vector w 
is w = V2 (eo + c-i + 6-2). The main difference is that now is the top 
eigenvalue of Fq with eigenvector v*. If w is chosen as in (|3.62[) . then 



(3.64) lira 2^F n s = (w | s) v* 

n »oo ^/g 



for all s e A4. The desired result 13.63|l follows when this is applied to s = 
/', •••/', e . □ 
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4. Wavelets on Fractals 

In section|21we proved some lemmas for representations of the Cuntz algebra 
O2, and we used them in section[3]abovc in our analysis of dyadic wavelets in the 
Hilbert space L 2 (R) . But there are other Hilbert spaces H which admit wavelet 
algorithms. In a recent paper |JJ, we constructed wavelets in separable Hilbert 
spaces H [s) built on Hausdorff measure (dx) s where s denotes the corresponding 
fractal dimension. For the middle third Cantor set, for example, the scaling 
number N is 3; and the fractal dimension is s = jr§ = log 3 (2). 

In this section, we show how our Cantor subdivision construction is built on 
a representation of the CVCuntz relations on L 2 (T) which are analogous to the 
representations of C 2 which we used in our analysis of dyadic L 2 (R)-wavelets. 

Lemma 4.1. Let N e N, N > 2, be given, and let m , m l7 • • ■ ,mjv-i be 
bounded measurable functions on T. Then the operators 

(4.1) (S i f)(z)=m i (z)f(z N ), feL 2 (T),i = 0,1,- •• ,N- I, z e T 

satisfy the Cuntz relations 



(4.2) 



S?Sj — S it jl 

N-l 

E s i s t = 1 



j=0 

if and only if the associated N by N matrix function 
(4.3) M N (z) = -±=(m j (ze i! &)) 



<N v v ' ' 

takes values in the unitary matrices for a.e. zeT. 

Proof. Rather than sketching the details, we will instead refer the reader 
to section |21 above where the argument is done in full for N = 2; see also [S] or 

m □ 

□ 

For each of the representations (Si) i= Q , we get a measure fio on the unit- 
interval [0,1], as outlined in section [21 (for the special case N — 2.) In the 
general case, the measure /j,q is determined uniquely on the iV-adic subintervals 
as follows, 



(4.4) ^ 



N N k ' N N k N k 



where i\, Z2, • • • € {0, 1, • • • , N — 1}. 

We state the next result just for the middle-third Cantor set; but, following 
the discussion in [2], the reader will convince him/herself that it carries over to 
any fractal constructed by iteration of finite families of (contractive) affine maps 
in R d ; so called iterated function systems (IFS). 
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Proposition 4.2. Let 
(4.5) < mi (z) = z 

Then the unitarity condition 14.3(1 from Lemma ETTl is satisfied, and the operators 



(4.6) (z) = m, (z) / (z 3 ) , / e i 2 (T) , i = 0, 1,2, z € T, 

define a representation of O3 on L 2 (T); and the corresponding measure Hq from 
(|4.4|) is the Hausdorff measure //W of Hausdorff dimension s = |^-| = log 3 (2) 
restricted to the middle-third Cantor set A3 c [0, 1]. 

Remark 4.3. A classical theorem of Hutchinson see also |10) . implies 
that the Cantor measure, supported on A3, is the unique Borel probability 
measure \x which solves 

(4.7) J j (x) dx (x) = ^ (/ / (|) ^ (x)+Jf ^±^j dfi (x) 

for all continuous bounded functions /. 

As a result, we need to verify that our measure /jq from 1(4.411 . and from the 
03-representation, satisfies identity 1)4. 7(1 . 

To compute the terms on the right-hand side in 14.4JI is a finite matrix 
problem. We need only to calculate the operators S*, i = 0,1, 2, on the three- 
dimensional subspace M.:— span {eo, e_i, e^}- But the argument from section 
[3] shows that the three restricted operators 



(4.8) F f . = S 

have the following matrix representation 



,z = 0,l,2 

M 




\ / 11 11 11 \ / 75 





^0 



















, F 2 = 




V 





/ 





(4.9) F = I -j= J,Fi=l 1 j,F 2 =l ^ 

/ V / V 

A substitution of 14. 8( 1- 14. 91 1 into 1(4.4(1 yields 



(4.10) ^0 



ii ifr ii if. 1 

N + " ' + A 1 ' N + '" + N T + N T 
\F ik ■ ■ ■F i2 F il e \\ 2 

if one or more of the ij 's is 1 
2~ k otherwise. 

Using this, it is immediate to see that fio satisfies Hutchinson's identity ((4.7(1 . 
and therefore //q is the Hausdorff measure // W , s = log 3 (2) , restricted to A3 . 
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5. A Technical Lemma 

In the proof of Lemma l3.t)l Theorem l3 . 1 2l and Corollarv l3 . 1 31 above . we relied 
on the following lemma regarding operators in a finite-dimensional Hilbert space. 
While it is analogous to the classical Perron-Frobenius theorem, our present 
result makes no mention of positivity. In fact, our matrix entries will typically 
be complex. 

Notation 5.1. If M is a complex Hilbert space, we denote by L(M) the 
algebra of all bounded linear operators on M If M is also finite-dimensional, 
we will pick suitable matrix representations for operators F:M — ► M If M 
contains subspaces, Mi, i — 1,2 such that 

Mi±M 2 and M = Mi ® M 2 , 

then we get a block-matrix representation 



(5.1) F = 



A B 
C D 



where the entries are linear operators specified as follows. 

A:Mi — B:M 2 — ► M x ; 

and 

C:Mi — >M 2 , D:M 2 — > M 2 . 
If dim .Mi = 1, and M\ = Cw for some w £ M, then we will identify the 
operators Mi — ► M with M via 

T V :C 3 z — ► zrj, 

where t) € M. The adjoint operator is 

T*x = {r) I x) w, for x € M. 

Lemma 5.2. Let M be a finite-dimensional complex Hilbert space, with 
d = dim M. Let F 6 L (M), and let a £ C satisfy the following four conditions; 

(i) a £ spec (F) ; 

(ii) \a\ > max{|A| | A £ spec (F) \ {a}} ; 

(iii) the algebraic multiplicity of a is one; and 

(iv) there is & w £ M, \\w\\ — 1, such that F*w = aw. 
Then there is a unique ^ £ M such that 

(5.2) (w | = 1 and F£ = a£. 
Moreover, 

(5.3) lim a - n F n x = (w | x) £ for all x £ M. 

n— >oo 

Lemma 5.3. There is a constant C independent of d = dimM and of x, 
such that 

(5.4) ||a- rl F"x -{w\x)£\\< Cn^ 1 max | - | s £ spec (F) \ {a}j . 
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Proof. [Lemma 5.2] Set 

M': = M Cw = {x G M \ (w \ x) = 0} . 

Then 

(5.5) M = Cw®M', 

and we get the resulting block-matrix representation of F, 

a 00 • • • 



(5.6) F- Xrf Q 

where a is the number in (i), the vector rj G M, and operator G E L (M. 1 ), are 
uniquely determined. As a result, we get the factorization 

(5.7) det(A-F) = (A-a)det(A-G) 

for the characteristic polynomial. Assumptions (ii) and (iii) imply 

(5.8) spec (F) \ {a} — spec (G) ; 

and in particular, we note that a is not in the spectrum of G. Hence the inverse 
(a — G) _1 is well defined, and (a — G) _1 € L (M') . We claim that the vector 

(5.9) ^w + ia-Gy 1 ^ 

satisfies the conditions in l|5.2|l . First note that (a — G) 1 rj G M' , so (w | £) = 

(w | w) = \\w\\ = 1. Moreover, 

F£ = aw + rj + G (a — G) 1 rj = aw + a (a — G) 1 r\ = <z£, 

which proves the second condition in l|5.2|l . Uniqueness of the vector £ in l|5.2() 
follows from l|5.8|l . Using the matrix representation l|5.6|l . we get 

a 2 00 • • • 
ar) + Oq G 2 



and by induction, 
(5.10) F n = 



a n 00 • • • 

a"' 1 ^ + a n - 2 Gri + ■ ■ ■ + G"- 1 !] G™ 

a™ 00 • • • 

(a"-G")(a-G) _ S/ G" 

Hence, if we show that 

(5.11) lim a- n G n = 0, 



• > oc 



then the desired conclusion (|5.3|l will follow. Using the matrix form (|5.10|) . the 
conclusion (15.31) reads 



(5.12) lim a- n F n = ( 1 a °° n "° 

n — >oo \ (a — G) 77 

In proving (|5.11(l . we will make use of the Jordan- form representation for G. 
Jordan's theorem applied to G yields three operators D,V, N £ L (M) with the 
following properties: 
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(1) D is a diagonal matrix with the numbers spec (F) \ {a} down the diagonal; 

(2) V is invertible; 

(3) N is nilpotent: If d - 1 = dim (M!) then N^ 1 = 0; 

(4) [N, D] = ND - DN = 0; 

(5) G = V(D + N)V~ 1 . 

Let x £ M.', and let n> d. Using (2)-(5), we get 

a - n G n x = Va~ n (D + N) n V- l x 



-2 



\— in 

i 



i=0 



Va^D^N'V^x. 



But the matrix a n D n % is diagonal with entries 

{a- n s n ~ l | s £ spec (F) \ {a}} , < i < d - 1. 
Using finally assumption (ii), we conclude that 

and the proof of (|B . 1 1|> is completed. □ 



(5.13) lim 



Proof. [Lemma 5.3] Let the conditions be as stated in the Remark. From 
the arguments in the proof of Lemma 15.21 we see that the two vectors on the 
left-hand side in (|5.4|l may be decomposed as follows: 

(5.14) a- n F n x = (w\x)w+(l- a~ n G n ) (a ~ G) _1 77 + ar n G n P M :x; 
and 

(5.15) (w\x}^= (w\x)w+{a-Gy 1 f]. 
Hence, the difference is in Ai', and 

\\ a - n F n x- (w I 

= a~ n G n (P M 'X- (a-G) -1 »?) 

( s n 1 

< Cn d - X max | - | s € spec (F) \ {a}j 
which is the desired conclusion 1)5. 4[l . □ 
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